A New Window Onto Quantum Chaos 
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In this article the statistical properties of symmetrical random matrices whose elements are drawn 
from a g-parameterized non-extensive statistics power-law distribution are investigated. In the limit 
as g — > 1 the well known Gaussian orthogonal ensemble (GOE) results are recovered. The relevant 
level spacing distribution is derived and one obtains a suitably generalized non-extensive Wigner 
distribution which depends on the value of the tunable non-extensivity parameter q. This non- 
extensive Wigner distribution can be seen to be a one-parameter level-spacing distribution that 
allows one to interpolate between chaotic and nearly integrable regimes. 
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Over the past two decades random matrix theory has 
proven to be a useful tool in the study of quantum chaos. 
Universal signatures of quantum chaos were found for 
systems whose corresponding classical Hamiltonians ex- 
hibit chaotic behavior, and they are reviewed in ||l|, ||, |j . 
The level spacing distributions derived from ensembles 
of random matrices ^, |5| have been shown to describe 
closely the statistics of fluctuations in the energy spectra 
of some chaotic systems |6|. This concept has also been 
extended to scattering ^\\, dissipative systems |^ , inter- 
series mixing in helium |10|| and in general the analysis 
of systems with mixed regular and chaotic phase space 
regions [ pT| . 

In the generic case of mixed phase space, the Hamil- 
tonian can be written as H = Ho + eHi, with Ho an 
integrable Hamiltonian and Hi a non-integrablc Hamil- 
tonian. Regular and chaotic regions coexist in the 
mixed phase space and no universal laws have been 
shown to hold. Therefore level spacing distributions 
between neighboring eigenvalues have been suggested 
1^, |l2|, |l3[ Q that are composed of the distributions de- 
scribing the regular Ho and chaotic Hi regions separately. 
The regular regions have classically integrable Hamiltoni- 
ans Ho whose eigenvalues are uniformly distributed and 
therefore the nearest-neighbor spacing statistics are de- 
scribed by Poisson distributions [|l2j. The eigenvalues 
of the non-integrable Hamiltonians Hi in the chaotic re- 
gions on the other hand are assumed to have an ergodic 
phase space region and are thus taken to be Gaussian dis- 
tributed. The nearest-neighbor statistics for the chaotic 
regions then follow the Wigner-surmised level spacing 
statistics. These composite distributions then contain 
one or more adjustable parameters that allow for an in- 
terpolation between the Poisson and the Wigner level 
statistics. 

The individual regions have Hamiltonians Ho, Hi 
whose eigenvalues are assumed to be independent. Thus 
an assumption of statistical independence of the random 



matrix elements (the eigenvalues) is inherent to the the- 
ory also as a consequence of a semi-classical approxima- 
tion. This assumption is admittedly invalid in even the 
simplest quantum case if tunneling between the dif- 
ferent regions is not neglected. 

In this letter the level spacing distribution of a generic 
H is obtained from statistically dependent random ma- 
trix elements corresponding to correlated eigenvalues. 
The importance of using non-extensive statistic is that 
the nearest-neighbor level spacing distribution allows one 
to interpolate from the nearly integrable to the chaotic 
regime as the non-extensive parameter q is varied. The 
Hamiltonian matrix of H to be considered is an x TV 
real symmetric matrix, and is invariant under orthogonal 
transformations. 

The extensive Gaussian random matrix theory can 
be generalized by examining the non-extensive, or q- 
parameterized entropy ^]. For systems with sta- 
tistically dependent (say, Hamiltonian matrix elements) 
variables the joint probability decomposition is 



P{H,,H,)^PiH,\H,)P{H,), 
which gives the pseudo-additive entropy 

Sq{Ht,Hj) = Sq{Hj) + Sq{Hi \ H j) 

+ {l-q)Sq{H,)Sq{H, I H,), 



(1) 
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It is known [|7| that the Tsallis entropy satisfies this 
condition, and the resulting probability will be of the 
power-law form. The g- logarithm is luqX = — ■^-j^^zjj- ■ 
In the limit as g — > 1 the usual form of the natural loga- 
rithm and thus the extensive statistics and its exponen- 
tial (Gaussian) distributions is recovered. 

The entropy to be maximized given the constraints is 
then 



{S)q = 



l-fPl,{H)dH 
(9-1) 
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{Trm)q = J Tr{H') P'^^{H)dH 
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(3) 

and which is subject to the extra normalization condition 



J PN{H)dH = 1. 



(4) 



Then maximizing (S)^ with the above constraints, 
yields the least biased probability density or non- 
extensive distribution function 

Pn{H) = An (1 + 0{q - l)Tr{H^))^ , (5) 
where is a normalization constant and 



Ht 



(6) 



l<i<Ar 



^<i<j<N 



This is then the Tsallis power-law form of the probabil- 
ity density for the N x N random matrix elements. This 
nonextensive distribution then is use to obtain the corre- 
lated distribution function for the eigenenergies assuming 
that the given Hamiltonian is symmetric, real and can be 
diagonalized by means of an orthogonal transformation. 
The correlated distribution is given as 
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Consider a special case of non-extensive ensembles of 2 x 
2 matrices. The nearest neighbor spacing distribution 
p{s) is obtained from the correlated energy distribution 
function P{Ei,E2) by 



dEi / dE2P{Ei,E2)6{s - \Ei - E2\) 

-OO J —OO 

dEi / dE2\Ei-E2\ 

-OO J —OO 

x(l + (3{q-l)J2Ef] S{s-\E,-E2\).{8) 



The constants A and /3 are fixed by the two normalization 
conditions / p{s)ds = 1, / sp{s)ds = 1. 

Integrating Eq.7. with the above conditions yields the 
non-extensive nearest neighbor spacing distribution. 



p{s) = 



(5 - 3g)/3 



s(l + /3(g-l)sV2)'-'" 



where 



IT 1 



■r(^-2) 



(9) 
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FIG. 2: P (s) Vs. s, q = 1.01. The non-extensive and ex- 
tensive Wigner distributions are nearly super-imposed. The 
Poisson distribution is plotted using a solid line. 



and 1 < 9 < |. 

The behavior of the non-cxtcnsivc nearest neighbor 
spacing distribution obtained by plotting its values for 
a range of the level spacing s, given the 'inverse vari- 
ance' (3 and the non-extensivity parameter q. In Fig.l. 
the q dependence of (3 is plotted for values of q be- 
tween 1 < q < 1.5. In Fig. 2. the Poisson (solid), ex- 
tensive (long-short dashing) and non-extensive Wigner 
(short dashing) distributions are plotted for a low non- 
extensivity parameter q value of g = 1.01. The non- 
extensive distribution is nearly superimposed on the ex- 
tensive Wigner distribution as is expected for q— > 1. In 
Fig. 3. the Poisson (solid), extensive (long-short dashing) 
and non-extensive Wigner (short dashing) are plotted for 
a high value of the non-extensivity parameter q = 1.38. 
Here the distribution is greatly shifted and approaches 
the Poisson level statistics distribution. 
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FIG. 3: -P (s) Vs. s, q = 1.38. The extensive distribution 
has been plotted with long-short dashing. The non-extensive 
Wigner distribution is plotted with the short dashes. The 
Poisson distribution is plotted using a solid line. 

In this letter Gaussian distributed N x N random Hamil- 
tonian matrix elements is generalized to the case of the 



non-extensive statistics and the resultant power-law dis- 
tributions. A derivation of the subsequent level spac- 
ing statistical distribution, the non-extensive Wigner dis- 
tribution, is given. This derivation is obtained from 
maximizing the non-extensive entropy of Tsallis of the 
N X N symmetrical random matrix elements. The re- 
sultant of nonextensive level-spacing distribution is q- 
parameterized and the nearest neighbor spacing distri- 
bution varies from near integrable to chaotic regimes as 
the non-extesive parameter q is varied. The major im- 
portance is the possibility of using non-extensive statis- 
tics to shed light on the quantum signatures of classically 
mixed systems. In future work it will be interesting to 
apply these results to Hamiltonians of mixed systems be- 
tween regular and chaotic regimes where deviations from 
the Wigner statistics become pronounced. 
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